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Abstract 

\ Let v(x, t) = v r e r + v eg + v z e z be a solution to the three-dimensional incom- 

pressible axially-symmetric Navier-Stokes equations. Denote by b = v r e r + v z e z 
Q^' the radial-axial vector field. Under a general scaling invariant condition on b, 

we prove that the quantity F = rv e is Holder continuous at r = 0, t = 0. As 
an application, we give a partial proof of a conjecture on Liouville property by 
Koch-Nadirashvili-Seregin-Sverak in [15] and Seregin-Sverak in [26]. As another 
application, we prove that if b G L°°([0,T], BMO~ l ), then v is regular. This pro- 
vides an answer to an open question raised by Koch and Tataru in [16] about the 
uniqueness and regularity of Navier-Stokes equations in the axially-symmetric 
case. 
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1 Introduction 

In this paper we study the three-dimensional incompressible axially-symmetric Navier- 
Stokes equations. In cylindrical coordinates, the velocity field v = v(x,t) is of the 
j_j ■ form 

v(x, t) = v r (r, z, t)e r + v e (r, z, t)e e + v z (r, z, t)e z . 



Here and throughout the paper, we write x = (xi,x 2 ,z), r = r(x) = \/ x\ + x\ and 

= e z (x) = 
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are the three orthogonal unit vectors along the radial, the angular, and the axial 
directions respectively. The radial, angular (or swirl) and axial components v r , v e and 
v z of the velocity field are governed by (see, for instance, |20j) 

' dtV r + b . Vv r _ (v^l + dr p = (A - ^)v r , 

(1.1) 

d t v 2 + 6 ■ Vt> 2 + <9 2 p = At> 2 , 
< b = v r e r + v z e z , V -b = d r v r + ^ + d z v z = 0. 

Here without loss of generality, we have set the viscosity constant to be unit. 

A special feature of the axially-symmetric Navier-Stokes equations is that the quan- 
tity T = rv e (x,t) satisfies an parabolic equation with singular drift terms: 

d t T + b • Vr + -d r T = AT. (1.2) 
r 

We remark that T enjoys the maximal principle. For this reason the axially-symmetric 
case appears more tractable than the full three-dimensional problem. 

Nevertheless, it is well-known that global regularity of the three-dimensional incom- 
pressible Navier-Stokes equations is still wide open even in the axially-symmetric case. 
But if the swirl component of the velocity field v 9 is trivial, independently, Ladyzhen- 
skaya [T7] and Uchoviskii & Yudovich [29] proved that weak solutions are regular for all 
time (see also [IB])- Recently, tremendous efforts and interesting progresses have been 
made on the regularity problem of the axially-symmetric Navier-Stokes equations with 
a general non-trivial swirl. For example, in (HE], Chen-Strain-Tsai-Yau proved, among 
other things, that the suitable weak solutions are smooth if the velocity field v satisfies 
r\v\ < C* < oo. Their method is based on the classical results by Nash [23], Moser [22] 
and De Giorgi [B]. In [T5], Koch-Nadirashvili-Seregin-Sverak proved the same result 
using Liouville type theorem for ancient solutions of Navier-Stokes equations. See also 
[26] for a local version. 

A velocity field is called an ancient solution if it exists in the time interval (— oo, 0], 
and it satisfies the Navier-Stokes equation in certain sense. A well known fact is that 
ancient solutions represent structures of singularity of evolution equations, which makes 
the study of ancient solutions an important topic. 

In this paper, we study the axially-symmetric Navier-Stokes equations under a more 
general assumption on the radial-axial velocity vector b. To be precise, we consider b 
such that 

b = b 1 + b 2 + b 3 , V ■ h = V ■ b 2 = V • 6 3 = 0, (1.3) 

where 

HSE(6 1 ) < b 2 = V x B, sup ||S|| B mo < C, sup r|6 3 | < C„. (1.4) 

-T<t<0 -T<t<0,x£M. 3 
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Some motivation and explanation for the condition and notations are in order. Here 
[— T, 0] is the time interval where a solution exists. We often take T = 1 for convenience. 
The number C* is an arbitrary positive constant and HSE(fei) is called "the hollowed 
scaled energy" , defined by 

HSE(6i)= sup Enih), E R (b 1 ) = sup ]- [ |6 X (-, t)\ 2 dx. (1.5) 

$<R<R -R?<t<0 ti Jb 2R /b r/8 

Here Rq is a positive number often taken as 1. 
We use 

||6|| E = HSE(6i)+ sup ||B|| BM o+ sup r|6 3 | (1.6) 

-T<t<0 -T<t<0 

to denote the controlling quantity of b throughout the paper. Here [— T, 0] is the time 
interval of concern, which may be shifted or scaled. The linear space consisted of those 
b such that < oo is called space E. The results in this paper depends on b only in 
terms of 1 1 & 1 1 js? - We will use a positive function .K"(||&||e) to denote such a dependence, 
whose precise value may change from line to line. Notice that the space E contains 
BMO~ x which is the largest known space in which the Navier-Stokes equations are 
well-posed. See the interesting work by Koch and Tataru [16J. Another feature is that 
the condition on b\ is imposed only on some subdomain of the space time cube. Outside 
of the subdomain, there is no restriction on b\. With a little bit more efforts, we can 
also just impose conditions on part of the space time for 6 2 and 63 too. But here we 
do not pursue that. 

Our first result states that T = rv e is Holder continuous at r = 0, t = if the 
radial-axial velocity field b satisfies (11 .3p and (jl.4p . The Holder continuity depends on 
b only through ||&| IE- 
Theorem 1.1. Given a number L > 0, let v = v(x, t), (x,t) G Ql = B(x, L) x [— L 2 , 0] 
be a L^ C (Q L ) weak solution to the three-dimensional axially- symmetric Navier-Stokes 
equations (jl.ip . Suppose that the radial-axial velocity field b satisfies (|1.3p -( J!T4"]) . Then 
T = rv e is Holder continuous at (0,0) uniformly, i.e. There exist positive constants a 
and C, depending only on ||6||e, such that, for all (x,t) G Ql/2, it holds 

\r(x,t) -r(o,o)| < c[(|x| + ^/\t\)/Lf sup|r|. 

Ql 

Our proof is inspired by [1] where the authors had proved a version of the above 
theorem under the assumption r|u| < C* using an De Giorgi type argument (see also 
[5] for the method based on the direct estimation of an evolution kernel). Here we will 
treat the more general b using a Nash type method in a uniform way. We will first 
establish a local maximum estimate for solutions of (11.21) in terms of the controlling 
constant C* for b in (jl.4p . This is done by using Moser's iteration method and De Giorgi 
type energy estimate, exploiting the structure of b. Similar argument has appeared in 
Zhang [30] and Chen-Strain-Tsai-Yau [H [5] where b is some form bounded function or 
\rb(x,t)\ < C*. Then we apply the Nash type method to prove the Holder continuity 
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of T. One handy tool which allows to treat more general type of vector fields b is a 
simple two dimensional integration by parts argument (12. 8p . Another tool is the John- 
Nirenberg inequality for BMO functions, which was first employed by Friedlander and 
Vicol |8j, and also by Seregin, Silvestre, Sverak, Zlatos |25j to treat the linear heat 
equation with Au + bVu - d t u = with b e L°°([0, T], BMO- 1 ). They prove Holder 
continuity of weak solutions to this equation. We also utilize the role played the stream 
function, which helps to do integration by parts one more time. Let v be a velocity 
field. We recall that a function B is called a stream function of v if v — V x B. 

The main significance of Theorem II. II is that it deduces the next two theorems. One 
of them gives a partial answer to an open question in [15] on Liouville properties. The 
other one establishes a condition on b such that solutions to axially-symmetric Navier- 
Stokes equations are regular. This regularity condition does not involve Lebesgue 
integral on b or absolute value of b, which may allow the capturing of more oscillatory 
functions. 

Theorem 1.2. Let v = v(x,t) be a bounded, weak ancient solution to ( II. ij) . Suppose 
also r\v e \ is bounded and the stream function is a BMO function. Then v = 0. 

Remark 1.3. The authors of [15] stated a conjecture on Liouville type theorem for 
the axially-symmetric Navier-Stokes equations: bounded, mild, ancient velocity fields 
are constants. The authors in [TJ)] proved such kind of Liouville theorems in the three- 
dimensional axially-symmetric case without swirl, or under the condition r\v\ being 
bounded. The above theorem, under the conditions that r\v e \ is bounded and the 
stream function is a BMO function, gives a proof of this conjecture. 

Recall that rv 9 is scaling invariant and it also satisfies the maximum principle. 
Therefore its boundedness is a mild restriction. A bounded function is obviously a 
BMO function. Although a bounded velocity field may not have a bounded stream 
function in general, a boundedness assumption on the stream function is also very 
mild since one expects it to hold in most natural cases when the velocity is bounded. 

In [IB] p25, Koch and Tataru wrote that there has been a strong interest in ob- 
taining well-posedness of Navier-Stokes equations assuming a BMO -1 space condi- 
tion. They raised the question of uniqueness and regularity for solutions in the space 
L°°([0, T], BMO -1 ). They also proved uniqueness for small solutions in such a space. 
Our Theorem 11.41 answers this question in the axially-symmetric case. We should men- 
tion that decay and analyticity estimates for solutions to the Navier-Stokes equations 
evolving from small initial data in BMO -1 are studied independently by Miura and 
Sawada [21] and Germain, Pavlovic and Staffilani [9]. 

Theorem 1.4. Let v = v(x,t) be a suitable weak solution to U.l\) in the space time 
region R 3 x [0, T\. Assume that the initial value satisfies v(-,0) E L 2 , \rv e (x,0)\ < C. 
Suppose also v(-,t) = V x B(-,t) with sup 0<i<r t) ||bmo < C*- Then v is smooth 
in R 3 x (0,T]. Here C and C* are arbitrary positive constants. 

Remark 1.5. Note condition \rv e (x, 0)| < C is only on the initial value. It can also be 
dropped by a approximation argument. We will not seek the full generality this time. 



4 



Remark 1.6. In [7J, Escauriaza, Seregin and Sverak proved that L^L 3 (Q) solutions 
to the Navier-Stokes equations are regular, which is the highly non-trivial borderline 
case of Serrin's criterion. Their proof is based on the method of backward uniqueness 
and unique continuation together with a blowup argument. Since L 3 is imbedded into 
BMO -1 , our Theorem 11.41 also provides a new and simpler proof to such a criterion in 
the axially-symmetric case. 

Before ending the introduction, let us mention some other related results on axially- 
symmetric Navier-Stokes equations. In the presence of swirl, there is the paper by J. 
Neustupa & M. Pokorny [23], proving the regularity of one component (either v r or 
v ) implies regularity of the other components of the solution. Also proving regularity 
is the work of Q. Jiu & Z. Xin [13] under an assumption of sufficiently small zero 
dimension scaled norms. We would also like to mention the regularity results of D. 
Chae & J. Lee [3] who prove regularity results assuming finiteness of another zero 
dimensional integral. On the other hand, G. Tian & Z. Xin [28] constructed a family 
of singular axially symmetric solutions with singular initial data; T. Hou & C. Li [12] 
found a special class of global smooth solutions. See also a recent extension: T. Hou, 
Z. Lei & C. Li [EE]. 

The paper is organized as follows: In section 2 we establish a local maximum 
estimate using De Giorgi type energy method and Moser's iteration method. Based on 
the local maximal estimate, we obtain the Holder continuity of T and prove Theorem 
II. II by Nash's method in section 3. The argument is based on [HE]. Then in section 4 
we prove our Theorem 11.21 and Theorem 11.41 using Theorem 11.11 and some new blow up 
arguments. The main idea is that a possible singularity falls only into two types. Type 

1 singularity can be scaled into an axially-symmetric, bounded, ancient mild solution. 
Type II can be scaled to a two dimensional ancient solution. Then we show that either 
type leads to a contradiction with the assumption that the stream function is in the 
BMO space. In the process the two dimensional Liouville theorem in [15J plays an 
important role. 

2 Local Maximum Estimate 

In this section we prove a local maximum estimate of T using Moser's iteration method 
in proving the parabolic Harnack's inequality. These estimates will be used to obtain 
Holder continuity of T in next section. The main idea is to exploit the divergence-free 
property of b(x, t) and to construct a special cut-off function. We also learned from 
[U [5] where the authors treated the term -d r T in the equation for T. 

We first derive an energy estimate of De Giorgi type for (11.21) . For this purpose we 
need a refined cut-off function. Set | < a 2 < (J\ < 1 and choose ipiy^s) = 4>(\y\)r](s) 
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to be a smooth cut-off function satisfying: 

supp0 C B(<Ji), (f) = 1 on B(a 2 ), < (j) < 1, 

supp^C (-(ai) 2 ,0], V (s) = lon(-(a 2 ) 2 ,0], < V < 1, (2.1) 
|„'| < i | v ^| < 1 IvC^I < l - 

. I ' I ~ (0-1-CT2) 2 ' Iv^l ~ 0-1-0-2' I \^/4>)\ ~ (01-0-2) 2 ' 

Here as usual we use A < B to denote the inequality A < CB for an absolute positive 
constant C . Such a cut-off function <fi can be simply chosen as a square of a standard 
cut-off function. We will also use the following notations for domains. Let R > 0, we 
write B R = B(0,R) and 

P(R) = B R x (-i? 2 , 0], P(R U R 2 ) = B Rl /B R2 x (-i? 2 , 0] for R 1 > R 2 . 

Consider the functions / = |r| 9 , q > | and the cut-off functions ip R (y,s) = 
My)Vn(s) = <j>(%)v(w)- Testin g O) by g|r| 2 ^ 2 r^| gives 



\JJ {dsf + (b ■ V)f + hrfyidyds = qJJ Ar\r\ 2 "- 2 r^ R dyds. 
Using Cauchy-Schwarz's inequality and integration by parts, we compute that 
q J j {/\T)\T\ 2 *- 2 Ti, 2 R dyds = qJJ (A|r|)|r| 2 ^V^ 

= ~ q II (( 2 ^- 1 )l vr l 2r29_2 ^ + v ^- |r| 29_1 v|r|)ti^s 

= - jj ((2 - i)|V/| 2 ^ 2 + 2^V^r • fVf)dyds 
= - JJ ((2-i)|V/|V> 2 +2fVij R - V(fij R )-2f 2 \V^ R \ 2 )dyds 
JJ \V(fij R )\ 2 dyds + JJ f\VM 2 dyds 



(2.2) 



< _ 



and 



l -JJ ^ R dJ 2 dyds 

= \f ^lf 2 (;t)dy -Iff f 2 d s ij 2 R dyds. 
Moreover, by the fact that T = on the axis r = 0, we have 

JJ -d r f 2 ip R dyds = JJ d r f 2 ip R drdzd9ds = JJ f 2 d r ip R drdzd9ds. 
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Consequently, using (12. ip . we have 

i>lf(;t)dy + If \V(fi. R ){'dyds (2.3) 
// ft* ~ \ //<» ■ ^ 



1 

2 

< 



Now we start to treat the drift term involving b = b\ + b 2 + &3- For i?i > R 2 , let us 
denote that 

E(R u R 2 ,b) = sup / |6(-,t)| 2 rfx. 

-i?2<i< -Ki - -K2 Jb Ri \b R2 

By ( 12. ip and the divergence- free properties of the velocity field 61 (x, i), we have 
~ // (&i • Vf 2 )^ R dyds =\l h- Y^ {M) l {r]R ff 2dyds 

L JJ J J P(<riR. ao.R) rhl 



P(<riR, cr 2 R) 

If 3 1 

E(aiR, a 2 R, &i)\5 



Therefore 

Next we treat the term involving 6 2 - Let B = B(t) be the average of B(-,t) in Br. 
Then 

-i // (6a • Vf)^ R dyds 

= 11 (B - B{t)) • V x {^^ R f)l( m ffz)dyds 

J JP(aiR. ao R) shl 



< 



+ II (B-B).[^xV^]dyds 



P(aiR, 02R) (j) R VVR 
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Therefore 



-\ff (62 " Vf)i, 2 R dyds 



< 



Hence 



L\L 2 (p(<jiK)) 

+ ((^)S II(B - B)/II ^(^)' 2 



-5// ft ' VfWldyds < i ^jT |V(,/>H/)| 2 A/<i 



To control the last expression, we need to recall the well-known John- Nirenb erg in- 
equality for BMO functions (see [H] or [27]): for any p G (0, 00), 

\\B(;t) - B(t)\\L H Bn) < C p \\B{.Mbmo\B r \ 1/p . (2.6) 
Taking p = 6 in the above inequality, we have 

^ll/il if ,3( P(CTlii) )PIU^M0|^|^i 

Plug this into (12. 5p . we deduce 



(2.5) 



-\ ff{b 2 -Vf 2 )^ 2 R dyds 

ll^ll Wog H/ii^ + W/ |V(«)| 2 ^ S . 

((Ji - a 2 ) 2 R 2 L*&{P(<nR)) 8 JJ P(aim 



PIIL-W* 8 ■■ 1 " (2 ' 7) 



The term involving 63 has been treated in [5] . Here we give an alternative proof for 
completeness and simplicity. 

\\ff (h ■ Wf 2 )^ 2 R dyds\ <\\JJ(h- Vmf)dyds\ 
= I JJ(h ■ ^^ R )M 2 )dyds\ < I J J ( l -\VM^Rf>drd6ds\ 



J I (\V^ R \^ R f 2 )drd6ds\ 

J I [d r {\V^R)f + \V^ R d r (f 2 )]rdrd9ds\. 



(2.8) 



Using Young's inequality, we deduce 

\\Jf{b,-Wf 2 )^ 2 R dyds\ 



^7 /„, \V^ R f)\ 2 dyds. (2.9) 



-a 2 ) 2 i? 2 " J H^(p{a lR )) ^ 8 JJ P{aiR) 

Plugging in the above three estimates (12. 4p . (12. 7p and (12. 9p on terms involving bi, 
i = 1, 2, 3 into (12. 3p . we arrive at 

sup [ ^ 2 R f 2 (;t)dy+ [[ \W{f^ R )\ 2 dyds (2.10) 

o-fR 2 <t<0 J B(aiR) JJPi^R) 



l + E{a l R ) a 2 R ) b 1 f ff r2 , , Wl^bmo^ 



~ (-i - °,) 2 R 2 JJpM f dyds + (-i - ^) 2 R 2 » J hi^ P{aiR) y 

By Holder inequality, this implies 

sup f ^ 2 R f 2 (;t)dy+ [[ \W{f^ R )\ 2 dyds (2.11) 

-a 2 R 2 <t<0J B(aiR) J J P{aiR) 



a 1 -a 2 ) 2 R^ U hfLs(p(v lR ))- 

Here and later in the section, as has been mentioned in the introduction, K = K(-) 
is a one variable function which may change from line to line, and is defined in 

Our next step is to derive a mean value inequality based on (12. lip using Moser's 
iteration method. By Holder inequality and Sobolev imbedding theorem, one has 
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(^ R f)*dyds< / f ||/^(-,s)||| 2(B( aiR)) \\V(MR)\\ 2 L 2 (B{aiR)) )ds 

P(viR) J ' 

4 

IIT7^„/,_M|2 

(o-ifl))- 



\ 2 (P 



-(o"i-R) 2 <s<0 

Using (EH]) and f l2TLj) . we obtain 







ff ,f dyds <! K »)\ ff fdyds \ 

J J P(u 2 R) (0"l — <y 2 ) Z R 2 J JPicrxR) > 

which implies that 

(in*)*** < { / (l|t|l ' )i p , // iri-4*}*. (2.12) 



P(a 2 R) L ((Tl - a 2 ) 3 R 2 J JP(cTiR) 
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For integers j > and a constant a = |, set er 2 = |(l + o"i = |(l + cr- 7 ); 

-- in fl2J2) . Then we have 

3(^y +i ^ e \ 3(lo) 




3 



yy f .(«(i+j.) 

By iteration, the above inequality gives 

, |rf^ +1 ^ 




3 



5 r-"'T>'- p // |r| 3 ^V 

it 2 J K JJp(R) J 



V 



We take the limit j — > oo to yield that 



sup |r| < (i^(||6|U)) 5 {^ / / iri 3 ^} 1 . (2.13) 



From this a well known algebraic trick (see p87 [10] e.g.) shows 



sup |r| < {K(\\b\\ E )){^ I f \T\ 2 dyds)\ (2.14) 



'P(R) 

Here the function K(-) may have changed at the last step. 



3 Holder Continuity of T 



In this section we study the regularity of T using the local maximum estimates of fl 2 . 1 3 [) 
in section 2 and Nash type method for parabolic equations. 

Let us first recall a Nash inequality, whose proof can be found in [5] . 

Lemma 3.1. Let M > 1 be a constant and \i be a 'probability measure. Then for all 
< / < M, there holds 



In / fd/j — / In fdjj 



< 



M\\g\ 



where g = In / — J In fdfi. 
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Let ( be a smooth radial cut-off function such that 

'C=l on B{\), ( = 0on5(lf, / R8 ( 2 (x)dx = 1, 



II<-> |v(?f)|< 



(3.1) 



oo, 



and C-r(^) = -\C(%)- Let $ be a positive solution to (1 1.2ft in P(-R). 

Lemma 3.2. Le£ $ < 2 6e a positive solution to (jl.2p m -P(-R) which is assumed to 
satisfy 

ll$ll ^(p(f))^ CojR5 " (3 - 2) 
Moreover, we assume that $(r = 0,z,t) a constant bigger than 1. Then there holds 

- y a(x)ln$(x ! t) C /x<M (l+||6|||) (3.3) 

/or a// i G [— £ x-,0] and some absolute positive constant M depending only on cq. 

Proof. First of all, let us define <3>(x,t) = &(Rx,R 2 t) and b(x,t) = Rb(Rx, RH). It is 
clear that $ solves the equation 

d t $ + b- V$ + -<9 r $ = A$ 

r 

on P(l) and < $ < 2, ||<&|| > Cq. The quantity we are going to control is 



— J CK^) l n = ~ j ( 2 ( x ) l n R~ 2 t)dx on a time internal [— £2 fs 0]. Equiv- 

alently, we just need to estimate — J ( 2 (x) ln $(x, t)dx for £ G [— y , 0]. 
Let \l/ = — ln $. It is easy to see that \l/ solves the equation 

d t ^ + b ■ W+ -<9 r # - A$ + IV^I 2 = 0. (3.4) 
r 

Hence, by testing (13.41) with ( 2 and using integrating by parts and Cauchy-Schwarz's 
inequality, one has 

d t J mC, 2 dx + J \V^\ 2 C 2 dx = J (-&• W--<9 r * + A^)C 2 dx 
<- / C 2 b-V(^ -y(s))dx - J ^d r (^ -^(s))C 2 rdrd6dz 
+ |Vv&| 2 C 2 dx + y |VC| 2 dx. 
Here ^(s) = / t)( 2 dx. Using the weighted Poincare inequality 

y |*-*(s)| 2 C 2 dz<<? y |V^| 2 C 2 dx (3.5) 
11 



and the divergence-free property of b, we can estimate 

1 



and 



C 2 6i • V(tf - V(s))dx 



(%■ V(V -V(s))dx 



< 



V*r( 2 dx + C / |6ir|VCrda; 



< 



< 



< 



VC 2 V x S(* - *(s))dx 
y |CV(*-*)||VC||S-S|dar + ^ |V(C^) ||S - B\\* - ^{s)\dx 

\(V(*-n\V(\\B-B\dx + J (t^ + y/fr^)\B-B\\C\\*-V{8)\dx 



\V*\ 2 ( 2 dx + C\\B\\l uo . 



Here we just used the weighted Poincare inequality and (12.61) . with p = 2. Moreover 

VC 2 ^ - ®(s))dx 
VC 



(%■ V(*-*(s))dx 
< /" |VC(*-*)|C^tZ0d2 



< | |CV(^-*)||VC|rfx+ J (^ + VCV^) |C| 
|Vv&| 2 C 2 rfx + C. 

Here we also used the integration by parts. On the other hand, by recalling the 
assumption that $(r = 0, z,t) is a non-zero constant, one can estimate 

—d r V( 2 rdrd6dz = -4tt / (* - ^)C 2 cfe + 4vr / (* - ^)d r ( 2 drd9dz 

r r=0 ./ 



47T 



^( 2 dz 



- 47T^> 



r=0 



C 2 ^ 



+ 4tt / O - mC—rdrdOdz 
r=o ./ r 



<C-C#(s) + ~ y \W\ 2 ( 2 dx. 



Here we also used the fact that the support of l\d r (\ is away from z- axis. Consequently, 
we obtain 

dtj v( 2 dx + c J ^( 2 dx<- 1 - J |w| 2 c 2 cfe + c(i + |HH). 

In order to proceed, we apply the Nash inequality in Lemma 13.11 Take / = $, 
d/j, = ( 2 (x)dx. One has 



In / dx 



mCdx 



$( 2 dx <M 2 \-m 



^( 2 dy\ 2 ( 2 dx. 
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Here M = 2 is the upper bound of $. Using the weighted Poincare inequality (I3.5p 
once again, we have 



In / <f>( 2 dx + / ^C 2 dx 



f <S>( 2 dx) 2 < C I |W| 2 C 2 dx. 



Hence, we finally obtain 



d^(t) + C ^(t)<C(l + \\b\\ 2 E ) 

- (2C)- 1 In I §( 2 dx + ^ 2 ( I $( 2 dx 



Let xi s ) be the characteristic function of the set 



W = { s6 [-I,0):||8 W || ll(Bj) >|}. 



By the assumption (13.21) and hence 1 1$|| L1 ( Dfl ^\ > Q), one has \W\ > ^f-. In fact, if 



N)) 



\W\ < *f, then 



^"^(Pd)) < / 2 l 5 (o)l^ 



11' 



c (2tt + l)co 
— as < < c , 



w 



c 2 



In / §>( 2 dx + ^ 



which contradicts with ( 13. 2ft . Thus, we have 

d&Qb) + C *(t) < C (l + ii&Hi) - scj-^W 



Note that the obvious consequence of this inequality gives 

t(s 2 ) < f + C e Co ( S2 - Sl ) (1 + ||6|||) 
for — \ < Si < s 2 < 0. Hence, if for some s G [—4, — — ) such that 



*(*„)<— (l+H&M +2 In 9 

Cq z 



then we are done since 



*(*)<*(«o) + ^(l + ||6|||) 



for all t G [sq,0). Otherwise, one has 



(l+IHU)+2|ln^ 
Cq i 



(3.6) 



(3.7) 
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for all s G [-|, -f). For s 6 If n [-|, -^), one has 

In / §( 2 dx > In / $dx > In — . 
J Jb 1 2 

Hence, by (I3.6p . we have 

Go 4 4 

Solving the above inequality gives 



(-1) 

The bound of the ^) depends only on c since f ) > and |W| > ^p. Starting 
from s = — ^ and using (13.71) . we have 

ns)<*{-j)+C e C »(l + \\b\\ 2 E ) 

for all s G [— ^F-, 0], which completes the proof of the lemma. □ 

As an corollary, Lemma [3.21 gives a lower bound of positive solutions of ( 11.2H . 

Corollary 3.3. Let Co and M 6e given in Lemma lcLM Then there exists a constant 
< 5 < 1 depending only on swc/j i/jcrf 

inf $ > -. (3.8) 

Kf) 2 

Proof. Using Lemma [3 .2\ we have 

M (l + ||6|||) >-J e R {x)\n<$>{t,x)dx 

= ~ Cr( x ) ki®{t,x)dx - / Cr(x) ln®(t,x)dx 

J&<$><1 J$<s 

~ / (r( x ) ln$(t,x)dx 
>0- / (l(x)]n<$>(t,x)dx-\n2 ( R (x)dx 

J<1><8 Jl<<&<2 

>-/ Cr( x ) ln$(t,x)dx - In 2, 
which implies that 

C|(x)ln$(t,x)da;<l+||&||| 

$<<5 
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for < t < 0. Consequently, we have 

\{xeB(^)mt,x)<8}\<^(i + \\b\\%) 

for — < t < 0. Using the mean value inequality (I2.14p . one has 



sup (<*-*)+ < I* 1 // («-*)l**} 4 



< — - — K 



which gives 



mf $>(5--^=lf(||6|| B ) 

for some Co > which is independent of S and R. Then (13. 8 p follows by choosing a 
sufficiently small 5 such that 

5<exp{-#(||6|U)}. (3.9) 

□ 

Now we are ready to give the proof of Theorem 11.11 
Proof of Theorem 11.11 

Without loss of generality we take L — 1. 
For < r < 1, we define 

m r = inf T, M r = sup T, J r = M r — m r . 

P W P(r) 



As in [5], we define 



' 2(Mi-r) 



if Mi > —mi, 



$ j^k) iiMl <-m l . 

It is clear that < $ < 2 is a non-negative solution of (H.2p in P(l) and a = $(r = 
0, z, t) is a constant bigger than 1. To verify that $ satisfies the condition (13. 2p . we 
need the following lemma on the lower bound of for < p < 1 as in [5]. 

Lemma 3.4. Suppose that b satisfies (ll.3p and (II .4p . Then for arbitrary p e (0, 1), $ 
defined above satisfies 

4||«f|L P(m|)) >C- 1 (#(||6|U))-"a. 
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Proof. Since the lemma is scaling invariant, we just take R = 1 in the proof. Let 
ip = 4>(\x\)rj(t), where G Cq° such that = 1 on Bi, <fi = on £>f, ^= and V^= are 
bounded, r] E such that 77 = 1 on [— |, — |] and r\ is supported in (—1,0). Let us 
test ([L2]) by P e (0, ~), to derive that 



y ^ (<9 S $ P + (6 • V)$ p + ^ r $ p )^ 2 ^s =P JJ A&Qp-^dyds. (3.10) 



Similarly as in [5], we have 

2 



r 



(d r <$> p )^dyds (3.11) 



> -C i I <5> p dyds+^a p . 



Here \y'\ = \/yJTyJ if y = (y h 8/2,1/3) ■ Likewise 

^ (-a^+pA^- 1 )^^ (3.12) 

2<F[^(<9^) + |VC| 2 - ^-^CAC]dyda- - yy |V($^)| 2 A/ds 
> -C yy _ 4 (P - !) yy |v($f^)|%rf s . 

Moreover, concerning the term involving b, we estimate it as follows: 

- y^ ^ 2 (&! • V)$ p dyds = JJ • V^ 2 dyds (3.13) 



> — C||frl||L?°L 2 fPl'l.i % lJ < £ P 



LfL2(P(l,i))ll* llLiL2(P(l,i)), 



and 



i)\b 2 ■ V)<S> p dyds = J J (B - B) • V x ($> p Vip 2 )dyds 
(B-B)-V x [($^) 2 ^]ci|/cis 



■0 



> -^II^^IL^^i^IIV^^IL^cpciI))-^ // \B-B\^dyds 



> 



R' 
2(p-l) 



yy iv($^)i 2 dyds-c yy 

- C( // $ 2 ^yrf S ) 1/2 ( f f (B — Bfdyds) 1 ' 2 . 
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By Holder inequality and ( 12. 6p . we have 
- j j ' 4,\b 2 -V)^dyds> 2 -^j^- J J |V($^)|%rf S -K(||6|U)||$|r L2p(p(l|)) (3.14) 



Just like (12. 9p . we also have 



tp 2 {b 3 ■ V)<S> p dyds > ^— ^ J j \V{<S>i^)\ 2 dyds -C Jj § p dyds. (3.15) 
Substituting fl3TT5j) . fl37Hj) . fl3~T3l . fl3~T2|) and fl3~TT|) into fl3~TUD . we deduce 



l^ p <C JJ ^dyds + CK{\\b\\ E )\\$T L 



L?L 2 (P(U)) 



<cK(\\b\\ E )\\n p L2P{P{hl) y 



which completes the proof of the lemma, since p G (0, 1/2) is arbitrary. □ 
Now we continue the proof of the theorem. 

By Lemma [3.4[ $ satisfies the assumptions in Lemma [3.21 for R = 1. By Corollary 
13.31 one h as 

inf $ > -. 

Noting that and m 1 < inf^^^ T < sup^^^ T < M 1; we have 

J^o = OSC/^r < (1 - ~) Jl (3.16) 
Iterating (13 . 16[) immediately shows that $ is Holder continuous at (0,0). □ 

4 Applications to axially symmetric Navier-Stokes 
equation 

This section is devoted to proving Theorems 11.21 and 11.41 We begin with 

Proof of Theorem 11.21 

By the assumptions of the theorem, we can apply Theorem 11.11 to deduce that the 
function T = rv e is Holder continuous at the space time point (0,0). More precisely, 
for any fixed point (x, t) e R 3 x (— oo,0), there exist positive constants a and C such 
that for all sufficiently large L > 0, we have 



|r(z,t)-r(o,o)| < c[{\x\ + v/|t|)/L] Q su P r. 



17 



Letting L — > oo, we find that 

rv $ (x,t) = r(x,t) = r(o,o). 

Since v e is a bounded function, the only way this can happen is v 9 = 0. Hence v is 
a bounded, weak ancient solution without swirl. According to Theorem 5.2 in [T5] . 
the ancient solution v = (0,0, l(t)) where I = l(t) depends only on time. Therefore its 
stream function B is a harmonic function since AB = -Vxd = 0. Since the function 
B = B(-,t) is BMO, by (EE]) we know 

/ \B(x,t) - B(t)\dx < C||-B('? t) \\bmoR 3 ■ 

J\x\<R 

Here Bit) is the average of B(-,t) in the ball Br. Since B(-,t) is harmonic, the mean 
value theorem tells us that B(t) = B(0,t). Hence 

/ \B(x, t)\dx < C\\B(-, t) \\bmoR 3 + 5(0, t)R 3 . 

J\x\<R 

The mean value theorem then implies that B(-,t) is a bounded function since 

\B(y,t)\ = -J— I B(z,t)dz 
4vr|y| 3 J B{yM 

/ \B(z,t)\dz<\\B(-,t)\\ BMO + \B(0,t)\. 

The classical Liouville theorem shows that the stream function B, being a bounded 
function, is constant. Therefore v — V x B — 0. □ 

Proof of Theorem 11.41 We use the method of contradiction. If there is a singularity 
to the axially-symmetric Navier-Stokes equations (11. ip . then we can generate a nonzero, 
bounded, weak ancient solution as in [TJ] . Our Theorem 11.11 and a scaling argument 
will then be used to show that such a bounded ancient solution is identically zero. This 
contradiction proves that singularity can not occur. 

By time shifting, we assume that the solution v exists in the time interval [—1,0] 
and that t = is a blow up time of v. The partial regularity theory in [2j [19] says that 
the Hausdorff measure of the singular space-time set of any suitable weak solution is 
zero. This implies that for axially-symmetric Navier-Stokes equations ( 11. ip . suitable 
weak solutions can only develop singularities on the symmetric axis r = 0. Hence, 
without loss of generality, we may assume that (0, 0) is the earliest blowup point. 

For k > 1, let (xk,tk) be a sequence of points such that 

-1 <t k /*0, Qk = \v(x k ,t k )\=j k max \v(x, t)\ /* oo, j k — > 1. (4.1) 
Define a sequence of functions {v^} by 

v (h) (x, t) = -Lv{x k + JL t fc + -Q 2 k (l + t k )<t< 0. (4.2) 
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By PP, one can assume that = r(xk) are uniformly bounded. It is clear that {v^} 
defined in f!4.2[) are mild solutions. Moreover, {v^} (up to a subsequence) converges 
to a bounded ancient weak solution u(x,t) to the Navier-Stokes equations (for details, 
see [15J). By the construction, \u(x, t)\ < 1 and \u(0, 0)| = 1. 

We consider two cases. 

Case 1 is when rk\v(xk, tk) \ = fkQk are uniformly bounded by some positive constant 
C. Then the functions {v^} are also axi-symmetric with respect to an axis which is 
parallel to the z-axis and is at distance at most C from it. Consequently, u is also 
axi-symmetric with respect to a suitable axis. Note that both the stream function 
and rv e are scaling invariant. Thus the stream function of u is in BMO and ru 9 is 
also bounded. Therefore we can apply Theorem 11.11 on u, which says that the swirl 
component of u vanishes. By Theorem 5.2 in [15], we conclude u = (0,0, /(t)) with 
I = l(t) being a function of time only. But this show u = as in the proof of the 
previous theorem. This contradiction shows that Case 1 can not happen. 

Case 2 is when Tk\v(xk,tk)\ = TkQk is not uniformly bounded. 

Hence, r k Qk (up to a subsequence) goes to infinity as k tends to infinity. Due 
to Caffarelli-Kohn-Nirenberg's partial regularity theory, {xk} (up to a subsequence) 
converges to which is a point on the z-axis such that r* = 0. Due to the axis 
symmetry of v, Xk can be chosen so that 9(xk) — > #oo for a O^. Hence, e r (xk) —> v and 
ee(xk) — > v 1 - = (— z/ 2 ,z/i,0) for an unit vector v = (vi, V2, 0). Here e r (x) and ee(x) are 
defined as in the introduction and (r(x),6(x)) is the polar coordinate of (x\,X2). 

It is clear that 

x k + ^ e B(x k , 7 §M=) for x E B(0,^Q k r k ), 
tk-( 7 ^ m ) 2 <tk + ^<t k <0 for -Q k r k <t<0. 

By the assumption on initial value and the maximum principle, we know 

\v%y)\<- for yeB(x k , r -±), t < 0, 
r k 2 

which shows 

l \^t)e e {x k + 4) I = l\v°{x k + ^t k + ^)\<^ ^ ^ 

( for (x,t) e B{0, y/r k Q k ) x {-r k Q k ,0}. 

Note that on B(0, y/QkTk) x (— QkTk-, 0], it is easy to see that e T {xk + -§-) — > v and 
eg(xk + 7T-) — > v 1 - as n — > 00. Moreover, for each k, is still a mild solution to the 

3D Navier-Stokes equations. By (]4.3}1 . there exists a subsequence of {v^} (we will still 
denote it by {v^}) and a bounded ancient solution u(x,t) to the 3D Navier-Stokes 
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equations on M 3 x (—00, 0], which is mild in the sense of [13], such that 
v {k \x,t) = -^-v r (x k + -^-,t k + Ja)er(xk + 77-) 
+ ^ (a* + 75-,** + + tH 

-> M = M V + uV + u z e z mL°°(n) 
for any compact subset H of f 3 x R and u(x, t) ■ v L = 0. Hence, 

u(x,t) = u r (x, t)v + u z (x, t)e z . (4.4) 
On the other hand, for (y, s) G B(x k , -7=^= ) x [*fc — ( ^q^ i *fc] > one nas 

-77- s ) e e(y) - s)er(j/)] 

= TT^t^d/) s)e r (2/) + AlA s)e e (y) + w 2 (?/, s)e«(y)] 

= d e [v^(Q k (y-x k ),Ql(s-t k ))] 

= Q k (d e y ■ V)«W (Q fc (j/ - x fc ). - **)) 

= Qfclx/I (e<>(2/) ■ V)v^ {Q k (y - x k ), Q 2 k (s - t k )), 



which gives that 



^r[-v r {x k + £-,t k + -^)e e {x k + (4.5) 

X u X 

+V 9 (x k + -pr-,t h + —^)e r (x k + —)] 
Wk W k V* 



Qfckfc + (e e {x k + ■ V)v {k \x,t) 



for (sc,t) G 5(0, y/Q k r k ) x (-Qj.r fc ,0]. Since r fc Q fc -» 00, we know Q k \x k + ^jH ->■ 00 
for fixed x. But the left hand side of (14.51) is bounded by definition of Q k . Hence, let 
k — > 00, we have 

(u x ■ V)u{x,t) = 0. (4.6) 

Note that the Navier-Stokes equations are invariant under rotation. Without loss of 
generality, we set v — e\ and v L = €2- Consequently, the limit function 

u(x, t) = u r (xi, z, t)e\ + u z (xi, z, t)e z , 
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is a bounded ancient solution to the 2D Navier-Stokes equations. By Theorem 5.1 in 
[15] , the limit 

u{x,t) = u r (t)v + u z (t)e z (4.7) 

depends only on t and that |u(0, 0)| = 1. By the argument in the proof of the previous 
theorem, the boundedness of the stream function of u in BMO norm implies that 
u = 0. This contradiction shows that Case 2 can not occur either. Therefore the 
assumption that v becomes singular at (0, 0) is false, proving Theorem 11.41 □ 
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